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ABSTRACT 

Measurements of plasma number density and electron temperature have been obtained in the plumes of 
laboratory arcjet thrusters using electrostatic probes of both spherical and cylindrical geometry. The two 
arcjet thrusters used in this investigation had different constrictor/nozzle geometries and operated on 
mixtures of nitrogen, hydrogen, and ammonia to simulate the decomposition products of hydrazine and 
ammonia. An increase in the measured electron number density was observed for both geometries with 
increasing arc power at a constant mass flow rate and with increasing mass flow rate at a constant arc 
current. For a given operating condition, the electron number density decreased exponentially off centerline 
and followed an inverse distance-squared relationship along the thrust axis. Typical measured electron 
temperatures ranged from 0.1 to 0.2 eV. 


INTRODUCTION 


The low power arcjet thruster is under development for auxiliary propulsion applications such as north- 
south stationkeeping onboard communications satellites.*^ As the arcjet thruster moves closer to being a 
flight ready propulsion system, spacecraft integration concerns become of significant interest. One key 
issue centers on the potential impact that the arcjet plume may have on the transmission of RF signals in 
common use for satellite communications. A first order study^ evaluated the interaction between a 1 kW 
arcjet plume and a 4 GHz communications signal in terms of signal attenuation and phase shift between 
transmitting and receiving antennas. In the antenna far field, the impacts appear negligible, except for 
propagation paths which pass very near the arcjet source. Numerical techniques are necessary to fully define 
the interaction processes and effects which would occur if the plasma were to expand within the antenna near 
field. Other integration concerns focusing on plume impacts include thermal loading and thrust losses due 
to impingement, particle contamination, and electromagnetic interference (EMI). 

A complete assessment of the potential plume impacts requires a detailed characterization of the arcjet 
exhaust. Initial research efforts^ demonstrated the applicability of Langmuir probes for the measurement of 
arcjet plasma properties and provided a partial description of the arcjet exhaust. Measurements 30 cm 
downstream of the thruster exit plane indicated plasma number densities on the order of 5 x 10^ cm‘3 at 
electron temperatures between 0.4 and 0.8 eV. The exhaust was only slightly ionized ( much less than one 
percent ). An increase in measured number density was observed with increasing arc current at a constant 
mass flow rate and also with an increasing mass flow rate at constant current. Radial measurements 18 cm 
downstream of the thruster exit plane indicated that the plasma number density decreased exponentially as a 
function of angle off centerline while values of electron temperature remained relatively constant with 
position. Similarly, an axial survey along plume centerline revealed plasma number density decreases as 
the inverse distance-squared. Finally, the initial study demonstrated an uncertainty in the true collection 
area of a spherical probe in a flowing plasma due to wake effects. 

This paper extends the previous work and presents a body of experimental data which correlate measured 
plasma properties in the arcjet plume with thruster operating conditions for two different constrictor/nozzle 
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geometries. The Langmuir probe is the primary diagnostic technique used in this investigation as it is a 
well established technique for the measurement of plasma properties. 

The first portion of this paper provides a discussion of the theory and operation of Langmuir probes for 
measuring plasma properties, along with the factors which govern interpretation of the probe characteristics. 
Details of the experimental hardware, data acquisition system, and test facilities are also provided. The 
major portion of this paper presents a body of experimental data on the plasma properties of the low power 
arcjet plume. Plasma number density and electron temperature are given as functions of arc power, 
propellant mass flow rate, mixture ratio, and constrictor geometry. A complete description of the plasma 
number density distribution is presented for a single operating condition. Measured values of the local flow 
angle with radial position are also discussed. Finally, a comparison of experimental and theoretical 
Langmuir probe characteristics and a study of the effective collection area of a spherical probe are given in 
the appendices. 


NOMENCLATURE 


Acff effective probe surface area, m^ 

Ac nozzle exit area, 
probe surface area, 
nozzle throat area, 
constrictor diameter, m 
nozzle exit diameter, m 
probe diameter, m^ 
electron charge, C 
arc gap setting, m 
current, A 

electron saturation current, A 
Boltzmann constant, J/K 
constrictor length, m 
cylindrical probe length, m 
charged specie mass, kg 
mass flow rate, kg/s 
number density, particles/m^ 

Maxwellian averaged cross section, cm 3/sec 
radial distance from center point of nozzle exit plane, m 
radius, m 


A P 

A* 

D c 

D c 

D p 

e 

G 

I 

loe 

k 

Lc 


ion speed ratio 


U, 


V 


2kf ~J~ 

/mj 


temperature, K 
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U G Mean drift velocity, m/s 

V voltage or potential, v 

z axial distance from thruster exit plane, m 
X dimensionless potendal (Eq. ( A1 ) ) 

X mean free path, m 

Xd Debye length, m 


Sub- or Superscripts 

A probe designation for fully conducting probe 

B probe designation for conducting/non-conducting sphere 

C probe designation for nonconducting/conducting sphere 

act actual 

e electron 

eff effective 

i ion 

max maximum 

min minimum 

n neutral particles 

p probe 

oo undisturbed plasma 


ANALYSES 


Langmuir Probe Theory and Operation 


The Langmuir probe is a widely used technique for the measurement of plasma properties. Briefly, it is 
a conductive element of known geometry which is connected to a power supply capable of biasing the probe 
to voltages both positive and negative with respect to the potential of the plasma to be investigated. The 
current drawn to the probe as a function of applied voltage yields fundamental information on the plasma 
properties, primarily electron number density and temperature. Excellent reviews of general probe theory 
and application in both the presence and absence of collisions may be found in Refs. 5 - 8. Laframboise^ 
established the theory for both cylindrical and spherical probes in a collisionless, stationary plasma. The 
application of the Langmuir probe to flowing plasmas ( both collisionless and collisional ) has been 
investigated by French Sonin* Clayden^, and Graf 13. 
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Provided that the probe docs not disturb the plasma ( or the macroscopic flow field ) and that the 
electrons exhibit a Maxwellian distribution, interpretation of the characteristic is reasonably straightforward. 
In the transition region, the electron current collected by the probe will be given by 1 1 


^e “ ^p e 


n z / \J 


kT c 

2jcm, 


exp 


<V P -Voo) 


kT„ 


( 1 ) 


Taking the logarithm of Eq. ( 1 ) and differentiating with respect to the probe voltage gives 


A in (I e ) = .-®_ 
dV e kT e 


( 2 ) 


Consequently, the slope of a semi-log plot of electron current as a function of probe voltage yields 
electron temperature. Plasma potential is generally located by extrapolating the linear portions of the 
transition and electron accelerating regions and finding the intersection. Using a thin sheath analysis, the 
saturation electron current ( at plasma potential ) is given by 


I 


oe 


= Ap 



(3) 


Rearranging Eq. ( 3 ) yields electron number density in terms of known quantities, as, 


n 


e 


(3.73 x 10 


'T ^oe 

Ap 


4 


11605 

T e 


(4) 


The above analysis is predicated on the assumption that the probe radius is small in comparison with 
the collision mean free paths but large in comparison with the Debye length; i.e. 


X » Rp » Xd 


where the Debye length is 


Xd 


= 6.93 x 10 
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(5) 


( 6 ) 


and X represents the collision mean paths. As long as the dimensions of the probe are small in comparison 
with the relevant mean free paths, the theory of Laframboise is a generally accepted means of determining 
the absolute currents collected by both spherical and cylindrical probes as functions of probe bias voltage, 
Rp/Xo, and Tj/T e . A method for predicting an ideal V-I characteristic in the collisionless regime is 
described in Appendix A along with a comparison of an experimental characteristic with an ideal 
calculation. 


Effects of Mass Motion 

In most plasma flows of interest, the electron random speed is much larger than the directed flow 
velocity. Previous investigators 1 1 1 >1 3 have found that cylindrical probes aligned with the flow 
streamlines behave as though they are in a stationary plasma. Spherical probes are promising because they 
do not have an alignment constraint; however, even when in free molecule flow with the neutral particles, 
spherical probe characteristics are influenced by the mass motion. The interaction is similar to that of an 
orbiting body in a rarefied ionized gas, 14 and is dependent upon potential, Ti/T e , and Rp/XD as well as 
ion speed ratio. Clayden 1 ^ has obtained reasonable results using spherical probes in a low density plasma 
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jet; however, a "wetted" current collection area ( one-half of the spherical surface area ) was used in the 
calculation of n e ( Eq. ( 4 ) ). This modification was made to account for the existence of a wake behind 
the probe which electrons were unable to penetrate, even though the probe was in near free molecule flow. 
A previous investigation of the arcjet exhaust using Langmuir probes 4 showed that there was some 
question of the true collection area of a spherical surface under varied flow conditions. Consequently, as 
will be described in Appendix B, a portion of the present effort was devoted to determining the effective 
collection area of spherical probes as a function of the ratio of probe radius to Debye length and the probe 
potential. 


Effects of Collisions 

Collisions between particles in a plasma influence the amount of current collected by a probe. In a 
flow field where the ionization fraction is small, collisions between charged and neutral species become 
important. For ion collection, the ion-ion and ion-neutral mean free paths are most important. Similarly, 
the mean free paths for electron-electron and electron -neutral collisions are relevant for electron collection. 
For the low density flows of interest here, it is not possible to maintain strict adherence to the above probe 
sizing criterion ( Eq. ( 5 ) ), particularly with regard to the ion mean free paths . Although the probes used 
in this investigation remained in free molecule flow with respect to neutral particles, they were not always 
small with respect to the charged species mean free paths. Table I lists the major collision mean free paths 
for the range of specie temperatures and densities encountered in this investigation. 

W ay mouth ' 5 has presented a theory for the operation of electrostatic probes under thin sheath 
conditions in plasmas in which the mean free paths are comparable with or smaller than the probe radius 
( i.e., RpAD > 100, X/Rp < 10, and X/Xq > 10 ). He found that the shape of the characteristic near 
plasma potential becomes distorted due to the probe perturbation of the plasma. The collisions reduce the 
current at the knee and at positive potentials and the reduction is sensitive to the ratio of Tj/T c . For thin 

sheaths, collisions in the sheath are negligible. The experimental data of Frenchl0 ( Sonin 11 , and Graf 1 ^ 
tend to support Waymouth's conclusions. However, all of the experimental data indicate that, as long as 
the probe is small compared to the electron mean free paths and that the flow speed is small compared to the 
electron thermal speed, the electron temperature may be found by the conventional means with a probe of 
any geometry, independent of its orientation. Further, French 1 ® has demonstrated that, under the conditions 
of both mass motion and small ion mean free path, the electron current may still be used to determine n e , 
as long as the mean free paths are much larger than the Debye length. 


EXPERIMENTAL APPARATUS AND PROCEDURE 


Arcjet Thruster Configurations 

Figure 1 displays a cross sectional schematic of the low power arcjet anode/nozzle configuration. The 
thruster discharge region incorporates a vortex stabilized, conventional constricted arc design. The tungsten 
insert shown in the figure serves both as anode and expansion nozzle. A complete description of the design 
features of the low power arcjet may be found in Ref. 1. The two laboratory arejets used in this 
investigation had internal engine configurations similar to that of the flight-type design but different 
constrictor dimensions. Table II lists the major components and dimensions of interest for the two arejets. 
The component materials for the two thrusters were identical except for the anode housings. For Thruster A 
( 0.038 cm D c ), the housing was fabricated from titaniated zirconiated molybdenum while, for Thruster B ( 
0.064 cm D c ), it was made from stainless steel. The different material selection enabled a larger operating 
range for the molybdenum thruster. Both of the anode inserts were fabricated from two-percent thoriated 
tungsten with a converging half angle of 30° and a diverging half angle of 20°. The arc gap setting for 
both thrusters was 0.058 cm. No attempt was made to vary this parameter for the present investigation. 
The area ratio of Thruster A was 600 while that of Thruster B was 220. 

Typical operating conditions and performance characteristics for the two arcjet thrusters may be found 
in Table III. Recent experiments 1 ^ have documented the relative effects of nozzle configuration on arcjet 
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operating characteristics. It was concluded that the contour of the divergent section of the nozzle/anode had 
a significant impact on arcjet performance, as it determined the location of arc attachment. A conical nozzle 
exhibited optimal performance. Since the anodes used in this investigation had the same divergent 
contour, the slight differences in measured performance were primarily due to plasma/gas dynamics in the 
constrictor region as it is not expected that arc attachment for the two configurations varied significantly. 
In general, Thruster B ( larger D c ) produced a higher thrust than Thruster A for the same ratio of arc power 
to mass flow rate although the overall efficiency and specific impulse were slightly lower. 

The arejets operated on mixtures of nitrogen and hydrogen to simulate fully decomposed hydrazine and 
ammonia. A 3:5:2 mixture of nitrogen, hydrogen, and ammonia was also used to simulate the gas mixture 
exiting a hydrazine catalyst bed. Typical mass flow rates ranged from 22 - 38 mg/s for Thruster A and 33 - 
45 mg/sec for Thruster B. The latter thruster was also operated with pure ammonia at 52-58 mg/sec. The 
gases were stored separately and mixed upstream of the thruster inlet. Digital flow meters, calibrated for the 
specific gas type, were used in conjunction with a flow controller to set and monitor individual volume 
How rates. Arc power was provided from a well-regulated ( 10 percent ripple ), pulse width modulated 
supply. The arc current was measured using a Hall effect current probe while arc voltage was measured 
with a high impedance divider network. The entire system operated with the arcjet anode tied to facility 
ground. Typical values of arc current ranged from 6 - 12 A, depending on the thruster configuration and 
mass flow rate. 


Langmuir Probes 

Langmuir probes of both spherical and cylindrical geometry were used in this investigation. The 
spherical probes were made of carbide steel bearing balls with radii of 0.951 and 0.554 cm. Each ball was 
silver soldered to a 0.76 mm tungsten wire. Stainless steel tubing was used to support the probe and 
alumina served as the insulator between the tube and probe. The cylindrical Langmuir probe was made of 
tungsten wire of radius 0.075 cm and length 2.38 cm and had the same type of support as the spherical 
probes. Further details of the design and construction of the Langmuir probes may be found in Ref. 4. 

Figure 2 displays an electrical schematic of the probe bias circuitry. The bias supply linearly ramped 
the probe voltage from a preselected negative value ( - 3 v ) to a positive value ( + 3 v) over a 2.2 sec 
time period. The probe measurements were referenced to the arcjet anode which was connected to facility 
ground. The probe voltage was measured directly across the output of the power supply while the probe 
current was obtained by measuring the voltage across a 101 Q resistor. Since the voltage drop across the 
resistor was not negligible, a correction was applied to the measured voltage during data reduction to obtain 
the actual probe voltage. The current and voltage signals were recorded by a digital storage oscilloscope 
and then sent to a personal computer for subsequent data reduction and storage. The above data acquisition 
process permitted a large number of complete probe characteristics to be taken during a given test period for 
later analysis. The estimated uncertainty in the probe current measurement was 50 - 100 jaA. 


Test Facility and Experimental Set-up 

The majority of the arcjet plume surveys were conducted in a 4.4 m diameter x 19 m long vacuum 
facility ( Tank 5 ). The vacuum tank was equipped with twenty, 32-in diameter oil diffusion pumps backed 
by four rotary blowers and four mechanical roughing pumps. It maintained background pressures less than 
2.6 x 10- 2 Pa ( 2.0 x 10' 4 Torr ) throughout the experiments. 

A photograph of the experimental apparatus and set-up is displayed in Figure 3. As shown, three 
Langmuir probes were mounted on an actuator system capable of linear motion in directions both parallel to 
and perpendicular to the thrust axis. The radial traversing mechanism covered distances of approximately 
± 0.5 m off the plume centerline while axial motion ranged from within 12 cm of the thruster exit plane to 
1.1m downstream. The cylindrical Langmuir probe ( center probe in the Figure ) was mounted on a rotary 
actuator such that it could rotate about a single point of reference for the purpose of studying local flow 
angle variation. The two spherical probes were mounted on either side and remained fixed with respect to 
the center probe's rotary motion. 
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Procedure 


The arejet thruster was mounted on a linear actuator which positioned the thruster inside the tank, as 
shown in Fig. 3. The probes were fixed at their predetermined locations prior to the experiment and the 
relative positions of the probes with respect to the center of the arejet nozzle exit were calibrated every time 
the vacuum tank was opened. A low power laser was used in conjunction with the center probe to ensure 
that the probe remained in the same vertical plane throughout its axial traverse. Axial and radial calibration 
of the actuators was accomplished with discrete measurements of the probe location, creating a linear plot of 
the actual versus indicated distance. Rotary motion of the center probe was calibrated in a similar fashion. 
The estimated uncertainty in an individual probe's position was ± 0.2 cm which is less than twenty percent 
of the spherical probe dimensions. 

The Langmuir probes were cleaned chemically before placement in the tank and every time the tank was 
cycled. It was not possible to heat the large probes sufficiently with a positive bias to get the temperature 
high enough for emissive cleaning. At least 2200 K for a few sec is needed to remove a monolayer of N2 

deposited on the surface of the probe. ^ An ohmeter was used to verify probe surface conductivity every 
time the tank was cycled. However, the effect of a changing and/or variable work function on a probe 
characteristic due to the presence of adsorbed N 2 on the probe surface will be discussed in a subsequent 
section of this report. 

Generally, the thruster was allowed to run for at least 30 min after the arc discharge was initiated so 
that it could reach a thermally stable condition. Subsequent changes to thruster operation, such as 
variations in arc current, required only a few minutes for the arejet to reach a new steady state operating 
condition. The arc current, arc voltage, propellant mass flow rates, and tank pressure were recorded 
periodically during the experiments. Previous measurements^ 30 cm downstream of the thruster indicated 
that, at a constant mass flow rate, n e increased significantly with arc power. For the present experiments, 
data were taken over a greater range of arc power and mass flow rates for mixture ratios of both 1:2 and 1:3 
of N 2 .H 2 . An axial distance of 30 cm was again chosen as the location for the measurements since the 
density levels at that location are sufficiently high to impact the transmission of RF communications 
signals. ^ 

Positional surveys were also taken at various locations in the plumes of Thrusters A and B. Radial 
surveys on both sides of the thrust centerline were taken at axial distances of 12.1 cm, 18.8 cm, 32.0 cm, 
and 64.5 cm. Axial surveys were taken along the thrust centerline and at the outermost radial location 
( 57.6 cm ). 

The current-voltage characteristic of a Langmuir probe was recorded by a digital storage oscilloscope 
and then transferred to a computer for storage and subsequent analysis. Approximately 2000 data pairs of 
voltage and total probe current were obtained for each probe sweep. Since the probe voltage was swept 
slowly, the total current collected by the probe reflected the temporal behavior of the plasma due to the AC 
ripple of the power supply; consequently, the total current data had to be smoothed with two, 16 point 
averaging routines. The fluctuations were generally less than 0.2 mA and were insignificant in the electron 
saturation portion of a V-I characteristic; however, they were the predominant cause of measurement 
uncertainty in the determination of electron temperature. Also, because the voltage was measured at the 
output of the power supply, the actual probe voltage was calculated by subtracting the smoothed current 
times the measurement resistance. The ion contribution was then subtracted from the total current to 
obtain the final set of probe voltage and electron current pairs. Least-squares fits were used to extrapolate 
the linear portions of the transition and electron accelerating regions of the characteristic. Plasma potential 
and the electron saturation current were defined at the point of intersection. The electron temperature was 
determined from the reciprocal of the slope of the line through the transition region and Eq ( 4 ) was used 
to determine n e . For the cylindrical probe data, the total surface area of the probe was used in Eq ( 4 ); 
however, for the spherical probe data, an effective collection area of 70 percent was used as will be 
discussed in Appendix B. 

In a flowing plasma, the magnitude of the saturation current depends upon the relative orientation of 
the flow and the cylindrical probe. If the ratio of Rp/^D »1 and the flow is supersonic, the ion saturation 
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current reaches a minimum when the probe is aligned with the flow streamline.^ A similar result is 
predicted for RpAD «1, as long as probe end effects are negligible.^ The dependence on probe 
orientation was utilized to measure the local velocity vector as a function of position in the exhaust. The 
probe was biased to + 3 V and the electron current to the cylindrical probe was recorded as a function of 
rotary angle with respect to the plume centerline. As will be shown in a subsequent section, the positive 
probe bias produced a minimum in probe current at the angle of the flow streamline. Data were taken on 
both sides of centerline at axial distances of 12.2, 16.0, and 18.5 cm. 


RESULTS AND DISCUSSION 

This section summarizes the results of the Langmuir probe surveys in the exhaust of a low power, dc 
arcjet thruster. The plasma plumes of the two different constrictor/nozzle geometries are described at a 
location 32 cm downstream of the exit as a function of arc power, propellant mass flow rate, and mixture 
ratio. Positional variations in electron number density and the local flow angle are also illustrated. A 
comparison of an experimental V-I characteristic of a spherical Langmuir probe with an ideal calculation is 
presented in Appendix A. Appendix B presents preliminary results in the study of the effective collection 
area of a spherical probe as a function of the ratio of probe radius to Debye length and the probe potential. 


Variations in Plasma Properties with Arcjet Operating Condition and Constrictor Geometry 

Figures 4 and 5 present cumulative plasma number density and electron temperature data for the two 
thruster configurations at a location on the plume centerline, approximately 32 cm downstream of the exit. 
The two thrusters operated on mixtures of N 2 , H 2 , and NH 3 to simulate the decomposition products of 
hydrazine and ammonia. The different symbols in the plots denote combined variations in mixture ratio and 
mass flow rate. Since the trends for the two thruster configurations are similar, they will be addressed 
together. Because a set arc current did not always give the a consistent arc voltage, the measured electron 
number density values are plotted as a function of arc power. The calculated uncertainties in the 1:2 and 1:3 
mixture ratio data mostly vary from 10 - 50 %, primarily due to the large uncertainties in relatively low 
values of measured electron temperatures. The data for pure ammonia with Thruster B varied in calculated 
uncertainties from 25 - 60%. However, within the limitations of the measurement technique, there is a 
reasonable consistency in the measured data. As demonstrated in a preliminary investigation,^ the measured 
electron number density increases dramatically with arc power at a constant mass flow rate. For Thruster B 
( Fig. 4 ( a ) ), the increase in n e with arc power appears smaller than for Thruster A ( Fig. 5(a)) 
although Thruster B had a much smaller operating envelope. There appears to be some noticeable scatter in 
the measured n e for Thruster A, particularly for the 1:2 mixture ratio at 34 mg/s; however, many of these 
data points were taken during the first day of testing when arc voltage fluctuations tended to be most severe. 
In general, a fifty percent increase in arc power results in a factor of 2-3 increase in the measured electron 
number density. For both 1:2 and 1:3 mixture ratios, there is only a slight increase in number density with 
increasing mass flow rate for a given arc power. The 1:3 mixture also produced higher values of n e than 
did the 1:2 mixtures, presumably due to the presence of more hydrogen. Finally, it is interesting to note 
that thruster operation with simulated, fully decomposed ammonia (1:3 mix ) resulted in much higher 
values of centerline electron number density than did the pure ammonia operation. 

Figures 4(b) and 5(b) display the corresponding values of measured electron temperature for 
Thrusters B and A, respectively. As in the preliminary investigation^, a systematic variation in electron 
temperature with arc current, arc voltage, or arc power is not evident, as the scatter in the data is fairly large. 
For the two thrusters tested in this investigation, the measured electron temperatures appear to be more 
consistent ( 0.1 - 0.2 eV ) but much lower than those of the arcjet tested previously ( 0.4 - 0.8 eV). It is 
apparent from the variations in and kT e with arc power and mass flow rate, that individual thrusters have 
distinct plume characters; however, the limited amount of data taken to date are insufficient to separate the 
roles that such variables as arc gap setting, arc anode attachment, constrictor/nozzle geometry, mixture ratio 
and strength of vortex propellant flow have on the charged species distribution in the exhaust. However, it 
is important to document the trends which do exist and establish a body of experimental data for subsequent 
analysis when the gas/plasma dynamics of the internal arcjet are better understood. 
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In general, the plasma plumes of the two arcjet thrusters exhibited the same experimental trends; that 
is, as long as the arcjet maintained a constant arc voltage for a given arc current, an increase in centerline 
electron number density was observed with increasing arc power at a constant mass flow rate or with 
increasing mass flow rate at a constant current. These trends are also evident in the data of Figure 6(a), 
which plots the measured electron number density as a function of specific power ( the ratio of arc power to 
mass flow rate ) for Thruster A. Some of the data taken during the first test date with this thruster have 
been excluded from the plot. As seen in the figure, an increase in specific power results in an increase in 
the measured electron number density at a constant mass flow rate; however, increasing the specific power 
at constant current (e.g., corresponding to a decrease in mass flow rate ) results in a decrease in n e . It is 
expected that the observed trends in measured plasma number density in the plume should be consistent 
with arcjet performance characteristics. Because the electron number density profile does not vary 
significantly with thruster operating conditions, the measured centerline electron number density is 
indicative of the relative amount of frozen flow losses. Figure 6(b) plots the overall efficiency as a 
function of specific power for Thruster A operating on a 1:2 mixture ratio of N2:H2 at two different mass 
flow rates. The observed decrease in overall efficiency with increasing arc current at a constant mass flow 
rate is qualitatively consistent with the plume electron number density trends. The trends for decreasing 
mass flow rate at constant current, however, are in direct contrast. It would seem that, at a higher mass 
flow rate, recombination would be enhanced by the higher pressure and the frozen flow losses would 
decrease. The performance characterization reflects that this may indeed be the case in that overall thruster 
efficiency does improve with mass flow rate at constant specific power. However, the electron number 
density increases under these same conditions ( Fig. 6(a)). One possible explanation 1 ^ may be that, at 
higher mass flow rates, the arc becomes more constricted; consequently, the value of centerline enthalpy 
may be higher resulting a higher centerline plasma number density. Further quantitative data near the exit 
of the arcjet nozzle are needed to fully the observed plasma plume anomalies. 


Plasma Number Density Distribution in the Arcjet Plume 

The spherical probes were used to obtain radial profiles at axial distances of 12.1, 18.8, 32.0, and 64.5 
cm downstream of the thruster exit plane. A centerline data point was also taken with the cylindrical probe 
at z = 12.1 cm. Axial surveys were also taken along the thrust centerline and at a radial distance of 57.6 
cm. Figure 7 displays a 3-dimensional plot of the measured electron number density distribution of 
Thruster B operating on a 1:2 mixture ratio of N2:H2 at a mass flow rate of 43.4 mg/s. The arc current 

was 10.2 A. Based on the measurements, the Debye lengths were on the order of 10" ^ to 10“ ^ cm so the 
ion density distribution should display a similar profile to maintain plasma neutrality. In a medium where 
the ionization fraction is very low ( much less than one percent ), as in the case of the arcjet exhaust, the 
ion-electron mean free paths are considerably larger than for charged specie-neutral collisions. Therefore, it 
is not expected that ambipolar diffusion rates determine the charged specie distribution in the plume. 
Rather, they are a direct consequence of the gas dynamic expansion of neutral species. The ions will tend to 
follow the neutral gas distribution due to high ion-neutral collision rates and hold the faster moving 
electrons back. 

Radial profiles were also taken with the thrusters operating on different mixture ratios, arc currents, 
and flow rates. Although not detailed in this report, the trends repeat the findings of the preliminary 
investigation.^ There is an exponential decrease in measured number density as a function of angle off 
centerline and an inverse distance-squared decrease in number density along centerline. 


Flow Angle Determinations Using Cylindrical Langmuir Probes 

The response of cylindrical Langmuir probes is extremely sensitive to the relative probe orientation 
with the flow streamlines. This dependence on probe orientation was utilized to measure the local velocity 
vector as a function of position in the exhaust. The probe was biased to + 3 V and the electron current to 
the cylindrical probe was recorded as a function of rotary angle with respect to the thrust centerline. Figure 
8(a) displays the measured electron current to the probe as a function of probe angle with respect to the 
plume centerline for the cylindrical Langmuir probe located on the thrust centerline, 18.5 cm downstream of 
the exit. The solid data points are from a survey taken 25 min later in time than the open symbol points. 
The minimum at 0° with respect to the thrust centerline denotes the local flow angle of 0°, as expected. 


9 



Values of local flow angle based on the cylindrical probe measurements arc plotted in Figure 8 ( b ) as a 
function of radial distance from centerline at axial distances of 12.2, 16.0, and 18.5 cm. The increase in 
flow angle with axial distance reflects the particle expansion. The same data are plotted in Figure 8 (c) as a 
function of angle off centerline. The solid line in the figure denotes radial flow. The fact that the measured 
flow angle and the geometric angle off the nozzle axis agree indicates that the flow is source-like in nature 
at dimensionless distances of z/D e > 10. 


CONCLUDING REMARKS 


Langmuir probes of cylindrical and spherical geometry have been used to investigate the plumes of low 
power, dc arejet thrusters operating on mixtures of nitrogen, hydrogen, and ammonia to simulate the 
decomposition products of hydrazine and ammonia. This paper extended previous work and presented a body 
of experimental data which correlate measured plasma properties in the arejet plume with thruster operating 
conditions for two different constrictor/nozzle geometries. It is apparent from the variations in n e and kT e 
with arc power, mass flow rate, and mixture ratio that arejet thruster plumes can be quite sensitive to design 
and operating condition. The data taken to date, however, are insufficient to determine the relative effects of 
such variables as arc gap setting, arc anode attachment, constrictor/nozzle geometry and strength of vortex 
propellant flow on the charged species distribution in the exhaust. The trends which do exist have been 
documented for subsequent analyses when the gas/plasma dynamics of the internal arejet are better 
understood. 

In general, the plasma plumes of the two arejet thruster configurations exhibited the same experimental 
trends; that is, as long as the arejet maintained a nonfluctuating voltage for a given arc current, an increase 
in centerline electron number density was observed with increasing arc power at a constant mass flow rate or 
with increasing mass flow rate at a constant current. A 1:3 mixture of N2:H2 produced higher values of n e 
than did 1:2 mixtures, presumably due to the presence of more hydrogen. Also, thruster operation with the 
1:3 mixture of N2-*H2 resulted in much higher values of centerline electron number density than did pure 
ammonia operation. Typically, the trends appeared to be more pronounced for the thruster with a smaller 
constrictor diameter. The measured electron temperatures consistently ranged between 0.1 - 0.2 eV for both 
thrusters but no systematic variation with arc current, arc voltage, or mass flow rate was observed. 

The measured electron number density decreases exponentially as a function of angle off centerline and 
as the inverse distance-squared along centerline. Radial profiles taken with the thrusters under different 
operating conditions produced similar trends. Finally, the use of a cylindrical probe for the determination of 
local flow angle was demonstrated. At dimensionless distances z/D e > 10, the measured local flow angle 
was equal to that of the geometric angle with respect to centerline, indicating that the flow was source-like 
in nature. 

It is apparent that more work is needed to fully document the arejet plume, particularly close to the 
nozzle exit. Spherical probes are promising because they do not have an alignment constraint but more 
work is needed to understand the presence of wake effects and the correct application of an effective probe 
area. Significant progress in the field of Langmuir probes has been made in understanding the idealized 
current collected by both cylindrical and spherical probes. Unfortunately, the theories are all based on an 
idealized, a priori knowledge of the major plasma parameters. It is much more difficult to extract 
information of the plasma properties ( such as ion temperature or ion mass) from an experimental V-I 
characteristic. A combination of intrusive probe and non-intrusive diagnostic techniques arc required to 
unambiguously obLain this information on the arejet plasma plume. 
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APPENDIX A 


Collisionless Theories for Electron and Ion Current Collection 


Theory and Development of Equations for Ideal Probe Characteristics 
In a collisionless plasma, the theory of Laframboise for both spherical and cylindrical probes was 
developed for a Maxwellian plasma at rest for the regime defined by 



0 < 

Tj/T e < 

1 


0 < 

RpAD ^ 

100 

and 

1 Xp Is 

25 



where Xp is the dimensionless probe potential given by 

Xp = e(V p -VJ/kT e (All 

The calculations of Laframboise assume that the dominant influence on the motion of the charged 
particles is the electric field consistent with Poisson's equation. For a two temperature plasma consisting 
of two species of charged particles ( electrons and singly charged ions ), the Vaslov equations governing the 
particle dynamics for the two species are coupled by Poisson's equation and solved iteratively. The absolute 
magnitudes of the ion and electron current densities collected by the probe as a function of probe potential 
may be summarized by the following relationships: 1 1 
For Xp < 0, 


For xp > 0, 


/kT e 
J i = e n^y — 

/ kT e 
Je -e noo y 2jtnic 

Ji = C n ~Y 

Je=e n ~V 2^; 


Ij (XpRp^'DT/I’e) 


P 



^e (XpRjADT/Fg) 


( A2) 
(A3) 


(A4) 


(A5) 


For Xp < -10, only ions are collected by a probe. For %p > - 10, both ions and electrons are collected 
by the probe. However, a widely used means of extracting the electron current from this portion of the 
characteristic is to extrapolate the ion current curve from Xp < ‘10 to higher potentials.^ In a plasma where 
the distribution of electrons is Maxwellian, the electron current collection in the transition region will 
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increase exponentially. For Xp > -1, the ion current collection becomes negligible until only electrons are 
collected. 

I*j and I* c in the above equations arc dimensionless ion and electron currents, respectively, and arc of 
order unity for moderate potentials. Kiel^O’^l an( j Peterson -Talbot^ have presented simplified methods for 
the calculation of I*i and I* e as functions of dimensionless potential, RpAD> and Ti/T e . In particular, 
Kiel^O developed an analytical expression for the saturation current regimes of a free-molecular spherical 
probe with a finite sheath. The expression is valid for Tj/T e > 0.01 and probe potentials large compared to 
the floating potential. Debye lengths must be sufficiently small to ensure formation of a sheath 

(R p Ad > 3 ). 

The current drawn to the probe as a function of voltage is dependent upon parameters involving the 
sheath thickness and the electrostatic potential at the edge of the sheath. To determine the sheath 
thickness, Kiel derived an empirical expression to fit the more exact numerical solutions of the potential 
fields. The electrostatic potential was obtained from a transcendental equation and depended only on the 
ratio of Ti/T e . For large probe potentials, Kiel's expression agreed reasonably well with the numerical 
solution of Laframboise. Complete details of Kiel's technique and the derivation of equations may be found 
in Refs. 20 and 21. The method of Kiel was used to generate ideal characteristics for comparison with 
experimental data. Since the calculated Kiel values of Ij* and I e * approach unity for small probe potentials 
( i.e., as Vp approaches Voo ), this method was used beyond the stated range of validity so a complete 
characteristic could be obtained for purposes of comparison. 


Comparison of an Experimental V-I Characteristic with an Ideal Calculation 


Figure 9 displays a comparison of measured and calculated Langmuir spherical probe response in terms 
of current collection as a function of probe voltage. The measurement was taken with the spherical probe ( 
Rp = 0.554 cm ) at a radial location 35.3 cm off centerline and 18.8 cm downstream of the thruster exit. 
The calculated values are based on ratios of Ti/Te of 0.01 and 1.0. On a logarithmic plot ( Fig. 9(a)), 
the agreement between experiment and theory is excellent. The linear variation of electron current over at 
least two orders of magnitude is indicative of a Maxwellian distribution of electron velocities. The 
measured electron temperature based on the inverse of the slope of a least squares fit through the transition 
region is 1578 K ± 580 K. The measured electron number density is 5.5 x 10*tycm3 with a measurement 
uncertainty of 28%. No correction was made to the probe area as RpAD for this case was only 47. When 
plotted on a linear scale ( Fig. 9 ( b ) ), the data more clearly show a marked difference between experiment 
and theory in the functional dependency of electron current with probe voltage. Here the measured electron 
current is compared against calculated values for Tj/Te = 0.01, 0.1, and 0.2. The deviation is due in part to 
the breakdown of Kiel's theory for values of dimensionless potential Xp ^ 3. However, there is also a good 
reason to suspect that variations in work function of the probe surface contributed to the discrepancy. 
Firstly, the probe voltage was swept slowly enough (2.2 sec) that thermal variations in the probe surface 
could occur, and secondly, the presence of N 2 monolayers on the probe surface may have altered the surface 
work function. For this particular characteristic, the sheath was thin, and the relevant collision mean free 
paths were all larger than the probe radius. It is therefore unlikely that collisions played a role in the 
variation of the measured I e curve from the ideal. Finally, since electron current collection in the 
saturation region is independent of ion mass, it seems reasonable to assume that the ion temperature at this 
location is much less than the measured electron temperature ( Tj < 0.2 T e ). 

Figure 9(c) presents a comparison of experimental and ideal calculated values of current collected by 
the spherical probe at negative potentials. Total probe current is plotted rather than ion current although the 
electron current may be assumed to be negligible for values of dimensionless potential Xp less than -10. 
For the ideal calculation, the description of ion current collection in this region is highly dependent on an 
accurate knowledge of the ratio of ion temperature to electron temperature and the ion mass. Since these 
two parameters are not known with absolute certainty, the ideal curves were generated for two different 
temperature ratios ( Tj/T e of 0.1 and 0.2 ) and ion masses ( 3 x 10"^ and 4 x 10"^7 kg ). It should be 
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mentioned that these values of ion mass are considerably lower than that which would be calculated from an 
average molecular weight of the gas mixture ( assuming the ion mass is equal to the average neutral mass 
of 1.77 x 10' 26 kg ). Using the higher ion mass resulted in electron currents which are much too small. 
The fact that the lower values were needed to give good agreement between theory and experiment may 
indicate that the ionized species in the plume are mostly atomic hydrogen ions; however, given the fact that 
the collected ion currents are on the same order as the measurement uncertainty, no definitive conclusions 
can be drawn at this time. Spectroscopic data are needed to determine relative specie populations in the 
arc jet plume. 
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APPENDIX B - Study of Probe Effective Collection Area 


Experimental Apparatus and Procedure 

Because of the uncertainty in the true collection area of a spherical probe in a flowing plasma, a portion 
of the present effort was devoted to the determination of the true collection area of a spherical probe in a 
flowing plasma. Three probes of equivalent spherical geometry but varied composition were constructed. 
The diameter of each probe was 1.91 cm ( 0.750 in ). The first probe, designated as Probe A, was fabricated 
from a solid carbide steel bearing ball. Probes B and C were constructed of a solid carbide steel hemisphere 
bonded to a ceramic hemisphere with a high temperature ceramic adhesive. The conducting half of probe B 
and the nonconducting half of Probe C faced toward the arcjet thruster exit. The individual probe 
characteristics could be uniquely determined using the same data acquisition system previously described. 
Additionally, the probe voltage could be varied manually with measurement of the actual probe voltage 
displayed on a digital multimeter. The current collected by the probe was also measured across a 101 ohm 
resistance by a digital multimeter. 

The experiments were conducted in a 1.5 m diameter by 5 m long vacuum chamber ( Tank 8 ) equipped 
with four, 0.82 m diameter oil diffusion pumps. The background pressure during the experiment was 
maintained below 6.7 x 10" 2 Pa ( 5 x 10~ 4 Torr ). Figure 10 presents a photograph of the experimental set- 
up for the effective probe area studies. As shown, the three probes were mounted 17.8 cm on a single rake, 
which was attached to a linear actuator capable of motion across the jet ( i.e., perpendicular to the thrust 
centerline). The arcjet thruster was mounted on a cradle connected to a push-pull rod that extended the arcjet 
into the vacuum tank. Once in place, it rested on a fixed table to ensure that the thruster position was stable 
and repeatable from one test date to the next. The push-pull rod was also used to vary the axial distance 
between the thruster exit plane and the probes. With this configuration, axial distances were limited to a 
range from 7.6 cm to 38.9 cm. Placement accuracy was ±0.3 cm. The radial positions of the probes in the 
flow field were given by a calibrated potentiometer mounted on the actuator assembly. The linear actuator 
allowed repeatable positioning of the probes on centerline to an accuracy greater than 0.02 cm over a 35.6 
cm range of motion. However, during actual operation of the vacuum facility, the probe positioning 
accuracy was significantly impacted by the presence of tank vibrations which resulted in a oscillatory probe 
motion. The positioning accuracy during the measurements was estimated to be ±0.5 cm. 

Thruster A was used in this experiment, operating at an arc current of 1 1 A with a resultant arc voltage 
of 1 17 v. The mass flow rate was 34 mg/sec using a 1:2 mixture of N2:H2. Probe data were taken in the 
following manner for the single thruster operating condition. The fully conducting probe (Probe A ) was 
positioned on the centerline of the exhaust. The axial position was set and recorded. The voltage across 
the probe was varied from -5 v to + 5 v in 1 v increments with the total probe current recorded at each 
voltage. The data were taken at several axial positions ranging from 38.9 cm to 7.6 cm. It was not 
possible to obtain data at distances greater than 38.9 cm without having the thruster support structures 
interfere with the plume expansion. Probe B was then placed on centerline and current data for all eleven 
voltages were taken at the same axial locations as Probe A. Finally, the same procedure was repeated for 
Probe C. Approximately forty-five minutes elapsed during data acquisition for the three probes. During 
this time the arc voltage varied less than one percent, so it is not expected that fluctuations in arcjet 
operating condition significantly affected the probe current collection. Complete V-I characteristics of the 
three probes were recorded at z = 30.5 cm. 

Since the current collected by a Langmuir probe is directly proportional to the true collection area, an 
effective probe area may be calculated using values of currents collected by the three probes while at the 
same potential. The two methods which were used to calculate the effective current collecting area of Probe 
A are given below. In both methods it was assumed that wake effects, if present, would only disturb the 
back hemisphere; therefore, for a spherical Langmuir probe in which the effective collection area is the 
actual surface area of the probe ( 47cRp^ ), the amount of current collected should be twice that of Probe B. 
The ratio of the sum of the currents collected by Probes B and C to their total effective collection area 
should equal the current density of a fully conducting spherical probe. Edge effects between the conducting 
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and non-conducting portions of Probes B and C were neglected. Method 1 can be expressed by the 
following equation: 


Ib+ ! c _ 2i b 

A e ff 


Similarly, the current density of Probe A, based on an effective collection area, should also equal that 
of a fully conducting probe. Method 2 is given by 


JA = ^B 
A eff A act 


(B2) 


Equations ( B1 ) and (B2 ) do not distinguish between ion or electron current, as the total probe current was 
measured; however, a previous study^ has indicated that when I Vp | > 1.0 v, the probe is operating in a 
saturation regime ( i.e., collecting only ions or electrons ). 


Preliminary Results and Discussion 

The three complete V-I characteristics of the different probe configurations yielded good agreement in 
measured number density when the effective probe areas were considered. The measured electron 
temperatures agreed to within thirty percent ( the measurement uncertainty ). The electron number density 
as a function of axial distance was determined from the complete V-I characteristic at z = 30.5 cm and from 
inferred values based on the electron saturation currents collected by Probe A at the various locations. 

The ratio of the effective collection area to the actual probe area was calculated using the methods 
described by Eqs. ( B1 ) and ( B2 ). The two methods provided reasonably similar results, although 
Method 2 appears to produce a slightly higher effective area than Method 1. The cause could be attributed 
in part to the presence of a small amount of ceramic adhesive that contaminated a section a few tenths of a 
millimeter wide on the conducting hemisphere of Probe C ( which was used in the data analysis of Method 
1 ). For | Vp | > 3 v, the two calculation techniques agreed to within four percent. For | Vp | < 2 v, 
the differences varied as much as seven percent. It does not seem likely that a distortion of the sheath at the 
edges of Probes B and C contributed to the discrepancy, since the probes were all operated in a thin sheath 
mode. 

Figure 1 1 presents the ratio of the effective collection area of the fully conducting probe to its actual 
surface area as a function of the ratio of probe radius to Debye length for both positive and negative probe 
bias voltages. The values of the effective collection areas are based on the average values calculated by 
Methods 1 and 2. For positive bias voltages ( Fig. 11 ( a ) ), the effective probe area decreases from 75% at 
RpAD of 100 to 55% as RpAD approaches 600. A similar trend is observed for negative bias voltages, as 
seen in Fig. 11 ( b ). For RpAD > 200, the probe is definitely in continuum flow with respect to the 
ion motion. For supersonic plasma flows, a shock consisting only of ions and electrons is predicted for 
these conditions. 10 Even at the closest data point, z = 7.6 cm, the probe was only in transitional flow 
with respect to the neutral particle motion (X/Rp ~ 2 ). Because of the probe size and the limitations of 
the positioning apparatus, it was not possible to obtain effective area data for RpAD <100. However, the 
data of a previous study^ seems to indicate that, for RpAD« 100, it is possible for a spherical probe to 
have an effective area which approaches the total probe surface area. It should be noted that the Langmuir 
probes used in the plasma plume surveys in Tank 5 , had typical values of RpAD much less than 200. An 
effective probe area of 70% was used in all data reduction, as it seemed to be a reasonable average for the 
range of RpAD encountered. The positional variation in electron number density decreases almost two 
orders of magnitude within 50 cm; consequently, a constant effective area for all measurements would not 
introduce an unacceptable error in the number density profiles, as long as one accepts a measurement 
uncertainty of at least 30 percent. 
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TABLE I. - Estimates of Relevant Mean Free Paths for Langmuir Probe Operation 


Assumptions: n c ~ n, 

n e ~ 5 x 10"^ n n 

n e ranges from 9 x 10 2 to 1 x lO^/cm 3 
T e ranges from 0.1 - 0.3 eV 
Ti ~ 0.1 T e 



^min 

^•max 


Collision Type 

£cm} 

(cm) 

Notes 

ion - ion 

5 x 10' 2 

3 x 10 1 

Based on Spitzer's relation 1 lj 22 

ion-elcctron 

7 x 10 3 

4 x 10 5 

Based on Spitzer's relation 1 1 *22 

electron-electron 

3 x 10° 

2 x 10 2 

Based on Spitzer's relation 1 **22 

electron-ion 

3 x 10° 

2 x 10 2 

^e-i ~ ^e-e * * 

neutral-neutral 

1 x 10 1 

1 x 10 3 

Gas kinetic theory 

ion-neutral 

3 x 10° 

3 x 10 2 

i. 11 

A i-n A n-n 

electron-neutral 

2 x 10 1 

4 x 10 3 

Q - 5 x 10-8 cm^/sec 4 


TABLE II. - Major Thruster Components and Dimensions of Interest. 


Thruster A 


Thruster B 


Component Materials 


Anode Housing 

Titaniated Zirconiated 
Molybdenum 

Stainless Steel 

Anode/Nozzle 

2% Thoriated Tungsten 

2% Thoriated Tungsten 

Cathode 

2% Thoriated Tungsten 

2% Thoriated Tungsten 

Injection Disk 

Molybdenum 

Molybdenum 

istrictor Geometry 

Constrictor Diameter, Dc, cm 

0.038 

0.064 

Constrictor Length, L c , cm 

0.036 

0.025 

Arc Gap Setting, G, cm 

0.058 

0.058 

:zle Geometry 

Exit Diameter, cm 

0.932 

0.952 

Throat Diameter, cm 

0.038 

0.064 

Half Angle, deg 

20° 

20° 

Area Ratio 

600 

220 
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Table III. - Typical Arcjet Thruster Operating Conditions and Performance Characteristics 


Mixture Ratio of N2:H2 
Total Mass Flow Rate, mg/s 
Average Molecular Weight 
Inlet Pressure, N/m^ 

Arc Current, A 
Arc Voltage, v 
Arc Power, W 

Power/Mass Flow Rate, kJ/g 
Thrust, g 

Specific Impulse, sec 
Overall Efficiency, percent 
Cold Flow Isp, see 
Cold Flow Thrust, g 


Mixture Ratio of N2:H2 
Total Mass Flow Rate, mg/s 
Average Molecular Weight 
Inlet Pressure, N/m^ 

Arc Current, A 
Arc Voltage, v 
Arc Power, W 

Power/Mass Flow Rate; kJ/g 
Thrust, g 

Specific Impulse, sec 
Overall Efficiency, percent 
Cold Flow Isp, sec 
Cold Flow Thrust, g 


Thruster A f = 0.038 cm ) 


1:2 

1:2 

1:2 

1:2 

35.7 

35.7 

23.8 

23.8 

10.6 

10.6 

10.6 

10.6 

6.9x10 s 

8.9xl0 5 

6.0x10 s 

7.6xl0 5 

8.0 

11.0 

8.0 

11.0 

123.6 

118.7 

110.7 

106.6 

989 

1305 

885 

1172 

27.7 

36.6 

37.2 

49.2 

16.2 

18.1 

11.6 

12.9 

455 

506 

486 

541 

35.1 

33.1 

30.9 

28.2 

116 

116 

117 

117 

4.2 

4.2 

2.8 

2.8 


Thruster 11(0? = 0.064 cm ) 


1:2 

1:2 

1:2 

1:2 

49.7 

49.7 

40.7 

40.7 

10.6 

10.6 

10.6 

10.6 

3.93xl0 5 

4.23xl0 5 

3.41xl0 5 

3.69x10 s 

8.0 

11.0 

8.0 

11.0 

113.7 

105.2 

106.6 

98.6 

909 

1158 

852 

1084 

18.3 

23.3 

20.9 

26.6 

18.9 

21.0 

16.1 

17.8 

380 

422 

396 

436 

36.7 

35.8 

35.0 

33.6 

116 

116 

115 

115 

5.8 

5.8 

4.7 

4.7 
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ANODE HOUSING 



Figure 1. - Cross sectional schematic ot arcjet thruster constrictor region. 


FACILITY 



Figure 2. - Electrical schematic of Langmuir probe drive circuit and data 
acquisition. 
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Figure 3. - Experimental configuration for plasma plume measurements in 
tank 5. 
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ELECTRON TEMPERATURE kTe, eV ELECTRON NUMBER DENSITY, cm-3 


m, 

mg/s PROPELLANT 



(a) Electron number density as a function of arc 
power. 



ARC POWER, W 

(b) Electron temperature as a function of arc power. 

Figure 4. - Cumulative data for Thruster B (Z = 32 
cm; plume centerline). 
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(a) Electron number density as a function of 
arc power. 



(b) Electron temperature as a function of arc 
power. 

Figure 5. - Cumulative data for Thruster A (Z » 
32 cm; plume centerline). 
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(a) Electron number density as a function of 
power/mass flow rate. 



(b)Thruster A performance data using 1 :2 
mixture. 

Figure 6. - Plasma plume and performance 
characteristics of Thruster A (mixture ratio 
of N2:H2 - 1:2). 
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Figure 7. - Electron number density distribution in an arcjet thruster 
plume. Thruster B, I - 10 A, m ■ 43.4 mg/s (1:2 mix). 
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MEASURED FLOW ANGLE, deg MEASURED FLOW ANGLE, deg MEASURED PROBE RESPONSE 

(Ip x 101ft), V 



-20 -10 0 10 20 
PROBE ANGLE WITH RESPECT TO PLUME 


CENTERLINE, deg 

(a) Measured probe voltage as a function of 
probe angle with respect to plume centerline. 
R « 0 cm. 



RADIUS FROM CENTERLINE, cm 
(b) Flow angle as a function of radial distance. 
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(c) Measured flow angle as a function of 
angle off centerline. 

Figure 8. - Determination of flow angle using 
Langmuir probes. 


DIMENSIONLESS PROBE POTENTIAL 

(c) Total probe current collection as a function of 
of dimensionless potential. 

Figure 9. - Comparison of measured and calculated 
Langmuir probe response. 
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Figure 10. - Experimental configuration for effective probe area studies. 
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